On the Applicability of the Ergodicity Hypothesis to Mesoscopic Fluctuations 
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We evaluate a typical value of higher order cumulants (irreducible moments) of conductance 
fluctuations that could be extracted from magneto-conductance measurements in a single sample 
when an external magnetic field is swept over an interval Bq. We find that the n-th cumulant has 
a sample-dependent random part ±{{<? 2 )}' 1 '' 2 \/a„B c / Bo, where ({g 2 )) is the variance of conductance 
fluctuations, B c is a correlation field, and a„ ~ n\. This means that an apparent deviation of the 
conductance distribution from a Gaussian shape, manifested by non-vanishing higher cumulants, 
can be a spurious result of correlations of conductances at different values of the magnetic field. 



The ergodicity hypothesis (EH) plays a crucial role in 
mesoscopic physics. It identifies mesoscopic 'sample-to- 
sample' fluctuations of a certain quantity with an appar- 
ently random but reproducible spread in the magnitude 
of this quantity measured in a single sample as a function 
of some controllable parameter, e.g., an external mag- 
netic field or a gate voltage. The EH is based on the 
idea that, since the fluctuations result from the electron 
interference in scattering from impurities^S^ changing 
the interference conditions by, e.g., sweeping the field 
over a wide range, one obtains a data set (usually called 
a magnetofmgerprint) statistically equivalent to that for 
sample-to-sample fluctuations. 

Since most of experiments are performed, by necessity, 
on a few samples while practically any theoretical ap- 
proach uses some form of the ensemble averaging over 
disorder, 4 the EH provides the framework for comparing 
experiment with theory in mesoscopics. Such a frame- 
work appears to work so well that it is now routinely used 
by default, in particular in experiments^LS aiming to 
determine the shape of the conductance distribution. 

Here we show that for the ergodicity hypothesis to be 
valid for analyzing an apparent deviation of the distri- 
bution from the normal shape in systems with large con- 
ductance, one needs the range of magnetic fields so large 
that it may not be reachable in modern experiments. We 
show that a non-Gaussian shape of the distribution ex- 
tracted from a single-sample magnetofmgerprint can be 
an entirely spurious result of the non-ergodicity due to 
residual correlations of conductances at different values 
of the magnetic field. For brevity, we consider only con- 
ductance distributions, although our method and results 
can be extended to any mesoscopic observable. 

Let us start with a brief description of theoretical pre- 
dictions on the conductance distribution. Any distribu- 
tion can be described by its moments or cumulants (irre- 
ducible moments). 9 Within the standard one-parameter 
scaling approach to disordered conductors^ the cumu- 
lants of conductance fluctuations in the metallic regime 
at temperature T = (i.e. for a fully phase-coherent con- 
ductor) are given b y 11 ) 12 ' 13 

{{g-))=A nd {g) 2 - n , (1) 



where A n d are numerical coefficients dependent on the 
cumulant order n and dimensionality d. Here (g) is the 
ensemble-averaged conductance (i.e. conductance aver- 
aged over all the realizations of disorder) measured in the 
units of e 2 /h. In a good metal (g) 7$> 1, and the higher 
(n > 2) cumulants in Eq. are small 14 in compari- 
son with the universal variance (the second cumulant), 

((g 2 )) = (g 2 ) — (g) 2 = A 2 d ~ 1. A characteristic fea- 
ture of the Gaussian (normal) distribution is that all its 
high cumulants vanish. Therefore, Eq. means that 
the distribution of conductance fluctuations is expected 
to be almost Gaussian. An experimental detection of de- 
viations from the Gaussian shape, characterized by the 
higher cumulants in Eq. represents a challenging task. 

This is even more so, when one realizes that at T ^ 
Eq. is only valid for sample sizes smaller that phase- 
breaking length L$. Thus for a quasi-lD wire with length 
L > Lj,, conductance fluctuations can be estimated as 
resulting from independent contributions of iV = L/L^ 
conductors connected in scries, which lead, in accordance 
with the central limiting theorem^ to a further suppres- 
sion of higher cumulants *i£ 

i9 n )) L>H ~ i9 n )) L ^ ("f J ■ (2) 

However small are the deviations from the Gaussian 
statistics described by Eqs. QJ and (J2J, they would be in 
principle experimentally detectable, had one been able 
to accumulate sufficient statistical data from uncorre- 
cted different realizations of disorder - a task which is 
in principle possible for semiconductor nanostructures by 
applying a large number of annealing cycles^ Our main 
point is that the failure of the EH invalidates an appar- 
ently simpler alternative (and the only one known for 
metallic wires) based on extracting statistical data from 
a magnetofmgerprint g(B). 

To explain how this occurs, let us formalize a proce- 
dure of extracting the moments (g n ) B of conductance dis- 
tribution from g{B). It is done by sampling it into a 
histogram, so that (g n ) B are determined by the follow- 
ing averages^ (where a smoothed histogram fs{g) is the 
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distribution function) : 

(g n ) B = I g n {B)— = / g"( B )f B (g)— . (3) 
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The cumulants {(g n )) B are then found from these (re- 
ducible) moments using the standard relations 9 -, see 
Eq. jnjl below. 

It is well known^that for B , B' » B c values of g(B) 
are correlated for different B with the co-variance 

k(B-B') = ((g(B)g(B'))) = ((g 2 )) x K (^^) , (4) 

where K(x) decays slowly (as a power-law, see Eq. (fTHf) 
below) from K(0) = 1 to 0. Its form and the value of the 
correlation field B c (which corresponds to one magnetic 
flux quantum through the region of coherence) depend 
on the dimensionality, temperature regime, size and ge- 
ometry of the system. The ergodicity has formally been 
proved^ in the limit Bq/B c — > oo. For a finite Bq, one 
needs to evaluate to what extent the correlations (0J con- 
tribute to ((g n )) B - To estimate a typical value of ((g n )) B in 
a given sample, we calculate its disorder-averaged mean 
square, {{{g n )) B ). 

In the leading order in B c /Bq <C 1, we find 



B P . 



(((9 n )) B 2 ) = ((9 n )) 2 +a n -^((g 2 )} 



a„ = n< 



Bo 

[K(x)] n dx 



(5) 



This equation is the main result of the paper. It shows 
that the ergodic hypothesis works only for the variance 
itself, as (((g 2 )) B 2 ) = A((g 2 )) 2 with the coefficient of pro- 
portionality A being very close to 1 for B c /Bq <C 1. For 
n > 3, the second term in Eq. (J5J represents a systematic 
error in extracting a true value of the cumulant from the 
magnetofingerprint. No sampling procedure in the sta- 
tistical analysis of a magnetofingerprint would produce 
error bars smaller than this systematic error. In spite 
of the small parameter B c /B$, this error term is dom- 
inant for measurements in highly conducting wires due 
to a very fast decrease of the cumulants with n, Eq. (T|. 
further augmented for samples with L > L^, Eq. |J5j. 

Thus even a perfectly Gaussian sample-so-sample dis- 
tribution can generate a visibly non-Gaussian one ex- 
tracted from the magnetofingerprint of a single sample, 
owing to the contribution of conductance correlations in 
Eq. JSJl. Experimental values of cumulants extracted in 
this way are sound only if they comfortably exceed the 
systematic error, ±((<? 2 ))™/ 2 y/a n B c /Bo- This is mani- 
festly not the case for experimental results of Ref. Q • On 
the other hand, in multi-parametric conductance fluctu- 
ations experiments, such as in Refs. 5 6], where realiza- 
tions are effectively changed also by a sweeping the gate 
voltage (and thus changing the Fermi wavelength) over 
some interval Vq , such a systematic error may be reduced 



down to ±((g 2 )) n/2 ^(B C /B Q ) (V c /V ), where V c is a cor- 
relation value of the gate voltage. 

Before deriving Eq. JSJl, let us recall some basic rela- 
tions between moments and cumulants: 9 



(9 n ) 




(6) 



The summation here is carried out over all the partitions 
{mi} = {mi, . . . , m n } of n into the sum n = Y^k=i krrik, 
with TOfc being non- negative integers. Note that in a dia- 
grammatic approaches cumulants are 'simpler' than mo- 
ments: the n-th cumulant of the conductance distribu- 
tion is represented by the sum of all the diagrams made of 
n conductance loops connected upon the averaging over 
disorder, while the n-th moment is the sum of both con- 
nected and disconnected diagrams, the latter resulting in 
the contribution of the lower-order cumulants in Eq. JJjJ . 
For magnetofingerprints, when the moments rather than 
cumulants are naturally defined, Eq. @, one can use the 
inverse 9 - of Eq. © to express ((g n )) B via (<?")„: 



{mi} 



(-i)'(p-i)! n 



i 



(9% 



where p = mi + to 2 + . . . + to^. 

Then a straightforward route to Eq. (JjJJl would be to av- 
erage over disorder the square of ((g n )) B above. However, 
such a route would be very cumbersome due to a large 
number of cancellations of non-leading contributions. 

Instead, we will use a so called moment-generating 
function, 9 defined for any well-convergent distribution 

/(<?) by 



F(t) = / 

J — ; 



e gt f(g)dg. 



(7) 



The moments are given by coefficients of Taylor's expan- 
sion of F(t), while the cumulants are generated by 



lnF(t) = ]>>"))- => ((g n )) = d?lnF(t) 



t=0 



(8) 



We use this expression to obtain directly the disorder- 
averaged variance of ((g n )) B - It is more convenient to deal 
with central moments of the distribution, i.e. with mo- 
ments of Sg(B) = g(B) — go(B), for which the disorder- 
averaged (Sg(B)) = 0. Such a shift of the distribution 
center does not affect the cumulants. To simplify nota- 
tions, we will still write (g n ) B instead of (5g n (B)) B . 

Then the variance var [«<?™» B ] = (((g n )) 2 ) - (((g n )) B f 
of sample-to-sample fluctuations of ((<?"))„ is given by 



var 



t[({g n )) B ] =d?d?(\nF(t)x\nF(r)) 



t=o 

T=0 



, (9) 



where the moment-generating function J2J for the distri- 
bution is given by 



F(t) = 



B+B 



e g(B)t 



dB 
~Bo~ 



(10) 
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We keep small non-Gaussian corrections only in the l.h.s. 
of Eq. @, representing the disorder averaging in the 
r.h.s. as that over the normal distribution of g{B): 



(A[g(B)}) = - Vg(B)A[g(B)} 



S = \J 'iBdB l g{B)k- 1 g{B l ) , Z = J Vg{B) e~ s , 

where k _1 is the resolvent of the correlation function 10}. 
Equivalently, it can be expressed in terms of the Fourier 
transforms g u and of a random g(B) and the correla- 
tion function k(B — B'), with = g* and = «*: 



f dw |9a,| 2 re a 



(11) 



To calculate the logarithm of the moment-generating 
function we use the standard replica trick2i: 

IriF — lim 



1 



N^O N 

Using also Eq. we represent the variance (JHJ) as 



var[{(g n )) B ]= Hm 



X 



(n) 
NM 



n,m^o NM 



(12) 



where 



v( n ) _ fin fin 



t=0 

T=0 



(13) 



B+B N 

* J J l\ B ° L\ B * 

in) 

To calculate Xj^L we perform the integration over Dg 
first, turning it into the Gaussian integral over the 
Fourier components of g, Eq. 



1 



— \th,e + Tl^e k | 



= exp ^ y dcu 

After the inverse Fourier transform for k u , we arrive at 

X$l = d?d? f exp [t 2 V BB + t 2 V b , b , + trV BB ,} 
J B 



where 
V BB 

V B 'b> 

Vbb> 



N 
~2 

M 



N N 



k (°)+E E < B i- B 3) 

M M 

2 k (°)+E E K&k-K 



k—l m=k-\-l 
N M 



^2^2^ -B' k ) 



t=o 

T=0 



(14a) 
(14b) 
(14c) 



The derivatives with respect to t and r in the limit 

(n) 

t, t — > in the expression for X N ^j above allows one to 
keep only the terms proportional to t n T m in the Taylor 
expansion of the exponent, which leads to 



X 



(«) 

NM 



(n! 



E 

i 



-i/n—21 
V BB ' 



V 1 V 1 
V BB V B'B' 



(n-2l)\(l\f 



(15) 



The expression above is entirely equivalent to one that 
can be obtained directly from diagrammatic techniques. 
The replica trick, as always in perturbative calculations, 
is an exact tool that serves to eliminate automatically 
spurious diagrams, in the present case those that con- 
tribute to reducible moments but not to cumulants. 

The next step of the calculation consists in selecting 
contributions of the lowest power in B c /B from Eq. i|15[l . 
Each term in the sum contains a product of n correla- 
tion functions (@J, each being dependent on a difference 
between two -B's from the multiple integral (|13|) . The 
integration will produce a factor of {B C J Bo) m , where m 
is the number of different pairs of -B's in the product. 

To illustrate this point, let us first consider the term 
with I = in Eq. (|15fl . proportional to the n-th power 
of the sum of NM different k's in Eq. (|14c|) . The lead- 
ing contribution comes from the sum of the n-th power 
of each term in Eq. I|14c(l , as each contains only one dif- 
ference between different B's; as they are all equivalent, 
re-labelling reduces the I = contribution in Eq. (|T51) to 



X 



(nO) 
NM 



NM 



f K n {B-B>) = ^a,M 

J B -DO 



(16) 



=i fc=i 



with a n given in Eq. J5}, since the multiple integral i|13|) 
is reduced to a single integral over x — (B — B')/Bq, 
proportional to B c / Bq. All other contributions from the 
I = term are proportional to NM Pn,m (where Pn,m 
is regular in the replica limit), but they contain more 
integrals over powers of K, each adding an extra small 
factor of B c /Bq. 

Now we show that the terms with I ^ in Eq. I|15|) do 
not contribute in the leading order in B c /Bq. Schemati- 
cally, their contributions may be represented as 

H n N ~ 21 (aN + E^-i))' (aM + £ m( m-i))' , (17) 

where E's with appropriate indices stand for one of the 
sums in Eqs. i|14|) . The number of independent integra- 
tions over k's is not smaller than the number of dif- 
ferent E's: therefore, contributions from Eq. I|17|) that 
contain different E's will be small in powers of B c /Bq. 
On the other hand, those which contain only one E, like 

Z N ~M l (a 2 NM) 1 or (aA0" /2S M(M-i) , would vanish in the 
replica limit in Eq. I|12|). the latter term vanishing due to 
the fact that each E is proportional to (at least) the num- 
ber of the terms in it, indicated by its index. 

Therefore, Eq. i|16|) gives the only non-vanishing in 
the replica limit term proportional to the first power of 
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B c /Bo. Substituting it into Eq. i|12|) . we obtain the result 
of Eq. Let us discuss its quantitative implications. 
First we express K{x) in Eq. <@J a o 2 '?' 1 ^ 1 ? 

iir(ax) = [1 + .t 2 ]~ 7 , x = AB/B c , (18) 

for the three types of the conductors with different values 
of the exponent 7 and the correlation field B c : 

(a) 7 — 1/2 and B c ~ 4>o/[ w x nrin{L v , L s }] for a long 
wire of width w^minl-L^, L s } < L, where 0o = he/c 

(b) 7 = 1 and B c ~ y' -j^-ffto/L 2 for a dot of size L 

connected to the reservoirs via leads with N channels 
and mean level spacing A at temperature T > NA 

(c) 7 = 2 and B c ~ 0o /L 2 for a dot or a short wire (with 
L < mm{L v , L s }) at T = 0. 

For these cases, coefficients a n in Eq. (J5J are given by 

( Brin^)] 2 (a) 

a n = < " r (« - h)V* (b) 



The appropriate power of the variance in the r.h.s. of 
Eq. (0) represents an inevitable systematic error in ex- 



tracting the n cumulant from a single magnetofinger- 
print. This can be rather large. The widest range of 
magnetic fields achievable experimentally is limited to 
Bq < 20T by the necessity to use noiseless supercon- 
ducting magnets. For the lowest temperatures in exper- 
iments with metallic wires, T ~ 10 lOOmK, the cor- 
relations field B c ~ 10 4- lOOmT 5 - 6 7 , in agreement with 
the estimates for type (a) above. A meaningful devia- 
tion of, say, ((g 3 )) from the combined prediction of Eqs. 
JU and J2J) (not speaking of a more accurate estimate of 
the third cumulant^ that shows it to be of order 1/g 2 
rather than 1/g) which would manifest a breakdown of 
the one-parameter scaling, or inadequacy of the model of 
nonintcracting electrons in disordered potential should 
therefore comfortably exceed the value of 0.15((<? 2 )). So 
far, reliable experimental data for the higher cumulants 
in the metallic regime fall well below this limit. 
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